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We study the synchrotron radiation of gluons by fast quarks in strong magnetic field 

<N 

^ produced by colliding relativistic heavy-ions. We argue that due to high electric conductivity 

of plasma, time variation of the magnetic field is very slow and estimate its relaxation time. 



We calculate the energy loss due to synchrotron radiation of gluons by fast quarks. We 
find that the typical energy loss per unit length for a light quark at LHC is a few GeV 
per fm. This effect alone predicts quenching of jets with p± up to about 20 GeV. We also 
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I show that the spin-flip transition effect accompanying the synchrotron radiation leads to a 

^ strong polarization of quarks and leptons with respect to the direction of the magnetic fleld. 

Ch Observation of the lepton polarization may provide a direct evidence of existence of strong 
magnetic field in heavy-ion collisions. 

> 

•£} I. INTRODUCTION 

o 

cn 

It has been recently argued [1] that the magnetic field B produced by colliding relativistic 

heavy-ions can be as large as eB ~ m^/h at RHIC and eB ~ Ibm^/h at LHC [1, 2]. This is 

^ comparable to the Schwingcr critical value eBc = rr? jh for a quark of mass m. In such strong 

• ^ 

^ fields many interesting perturbative and non-perturbative phenomena are expected to take place 

H 

(see e.g. [1, 3-7]). In this letter we discuss the energy loss and polarization of fast light quarks 
moving in external magnetic field [8, 9]. In QED these phenomena have been studied in detail 
due to their significance for collider physics, see e.g. [10, 11]. In turn, energy loss of fast particles 
in heavy-ion collisions is one of the most important probes of the hot nuclear medium [12, 13]. 
Synchrotron radiation in chromo-magnetic fields was discussed in [14-16] . 

A typical diagram contributing to the synchrotron radiation, i.e. radiation in external magnetic 
field, by a quark is shown in Fig. 1. This diagram is proportional to (eS)", where n is the number of 
external field lines. In strong field, powers of eB must be summed up, which may be accomplished 
by exactly solving the Dirac equation for the relativistic fermion and then calculating the matrix 
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element for the transition q — )• q + g. Such calculation has been done in QED for some special 
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FIG. 1; A typical diagram contributing to the synchrotron radiation by a quark, 
cases including the homogeneous constant field and can be readily generalized for gluon radiation. 

II. TIME DEPENDENCE OF MAGNETIC FIELD 

At first we would like to determine whether the constant field approximation is applicable to 
magnetic fields created in relativistic heavy-ion collisions. For periodic fields, the measure of how 
fast the field varies is the Keldysh "adiabaticity" parameter 7 defined as [17, 18] 

where lob is frequency and is amplitude of the magnetic field. We can adopt 7 for adiabatically 
changing fields as well. In this case ujb = \B/B\ is the rate of the field change and Bm is average 
field over time I/lo. The constant field approximation corresponds to 7 ^ 1. 

Let us now estimate 7. At first, we assume that B is determined only by the valence charges 
of initial nuclei and neglect the magnetic response of the produced nuclear matter. In this case B 
decreases with time t according to the power law implying that ojb ~ 1/^ [1]) where t is the time 
measured in the center-of-mass frame. Obviously, at early times ujb is largest. Therefore, to set the 
upper bound on 7 we need to estimate the time to after a heavy-ion collision when the quarks are 
released from the nuclear wave functions. If Qs is the saturation momentum, this time is to ~ ^/Qs- 
For a semi-peripheral collisions of Gold nuclei at RHIC we obtain ujb = Qs ^ ^ GeV. Thus, the 
adiabaticity for u-quark is 7 = 0.1 — 0.2, while for d-quark 7 = 0.5 — 0.8 which marginally justifies 
the constant field approximation. Here we used = 1.5 — 3.3 MeV and = 3.5 — 6 MeV [19]. For 
heavier quarks 7^1 and it seems that we cannot apply the constant field approximation. Energy 
dependence of the magnetic field follows from its transformation properties under boosts eB oc , 
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where Y is the rapidity. On the other hand, Qg oc e^V^, with A < 1 [20, 21]. Therefore, we expect 
that 7 will decrease with energy improving the applicability of the constant field approximation. 

So far we have neglected the magnetic response of the quark-gluon plasma. The plasma seems 
to be strongly coupled and we expect that its dynamics in strong external magnetic field is highly 
non-trivial. Solution of the problem of plasma magnetic response requires extensive numerical 
simulations of the relativistic hydrodynamic equations. Bearing this in mind, we, however, can 
use semi-classical arguments to derive a simple estimate of how the time dependence is affected 
by the plasma magnetic response. We will denote the magnetic field due to the valence quarks as 
Bq. According to the Faraday law, decreasing magnetic field Bo induces electric field E circulating 
around the direction of magnetic field. The electric field generates electric current that in turn 
produces the magnetic field Bi pointing in the positive z direction according to the Lenz rule. In 

— * 

the adiabatic approximation the resulting total magnetic field B satisfies the following equations 
[22] 

V^B = nadtB, V-B = 0. (2) 

Here /x and a are electric permeability and conductivity of plasma. The initial condition sd, t = to 
reads 

B{f,to) = Bme-1^ , (3) 

where we used the cylindrical coordinates {z, p, (f)} and R is of the order of nuclear radius. We 
neglect the external magnetic field at t > to- 
Solution to the problem (2), (3) is 

B{r,t) = j dV'B{f^,to)G{r-7^,t-to). (4) 

where 

'^^'''')=(4;rtM372^^4"4t7-] 
the Green's function and we assumed /x = 1. Eqs. (4), (5) describe evolution of the initial field (3) 
in time. Integrating over the entire volume V' we derive 



B{r, t) = Bm „g , — im~ ^^P 



i?2 + 4(t-to)/a^L B? + A{t-to)/u_ 
It follows from (6) that as long as t — to ^t, where r is a characteristic time 



(6) 
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the magnetic field B is approximately time- independent. 



To estimate r we use the value of the electric conductivity found in the lattice calculations 
a ICemT'^/Tc [23], where Cem = ^TraEMY^f^}- We have a = 0A3Tc{T/Tc)^. For = 
170 MeV this becomes a 0.37 fm^^ [T/Tcf'. With a conservative estimate of the medium size 
i? ~ 5 fm we have r ~ 2.2 fm (T/Tc)^. Substituting ojb = 1/t into (1) we obtain 7 < 0.1 even 
for T = Tc- Moreover, recent lattice calculations show that the electric conductivity increases in 
the presence of strong magnetic field [24]. As plasma expands it cools down; at T < Tc electric 
conductivity rapidly decreases and magnetic field vanishes. We thereby conclude - keeping in mind 
the qualitative nature of our argument - that the magnetic field can be considered as approximately 
stationary during the process of synchrotron radiation. 

We could have concluded that the magnetic field is quasi-stationary by mere examining the 
diffusion equation (2). Indeed, on dimensional grounds magnetic field significantly varies over the 
time r ^ R^a, cp. (7). This argument is arguably qualitative. A quantitative analysis must include 
a more realistic contribution of external sources, effect of plasma expansion and realistic transverse 
geometry. The contribution of external sources would appear in (4) as an additional term involving 
integration over space and time (from to to t) of V x je, where je is the current density of valence 
quarks of receding nuclear remnants. Such corrections will certainly introduce time variations of 
the magnetic field. However, since r is more than an order of magnitude larger than the plasma 
life-time, we expect that adiabatic approximation 7 <C 1 still holds in a more accurate treatment. 
Of course, as has been already pointed out, a comprehensive numerical simulation is required in 
order to obtain the precise time and space dependance of the magnetic field. 

In light of our discussion in this section, we will treat magnetic field as slowly varying in time. 
Formulas for energy loss in Sec. Ill should be understood as time average, and S as a mean 
magnetic field. On the other hand, polarization effects discussed in Sec. IV, are sensitive only to 
the initial value of the magnetic field B{to). 

III. ENERGY LOSS DUE TO SYNCHROTRON RADIATION 

Quark propagating in the external magnetic field radiates gluons as depicted in Fig. 1. The 
intensity of the radiation can be expressed via the invariant parameter x defined as 





{pxBf 



(8) 
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where the initial quark 4-momentum is p'^ = (£,p), is the quark charge in units of the absolute 
value of the electron charge e. At high energies, 

Bcvn 

The regime of weak fields corresponds to % <C 1, while in strong fields x » 1. In our case, 
eB /eBc ~ (m^/m„)^ S> 1 (at RHIC) and therefore x 3> 1. In terms of Xi spectrum of radiated 
gluons of frequency oj can be written as [8] 

Ai«)^+(?.^,.'/=)Ai'w}. (10) 

where / is the intensity, 

and e' is the quark's energy in the final state. Ai is the Ayri function. Eq. (10) is valid under the 
assumption that the initial quark remains ultra-relativistic, which implies that the energy loss due 
to the synchrotron radiation Ae should be small compared to the quark energy itself Ae ^ e. 

Energy loss by a relativistic quark per unit length is given by [11] 

de p dl n?x^ p 4 + 5x^3/2 + 4x^x3 

— = — / dLO—=asCF / ; TTT^A Alixjxdx. (11) 

dl Jo dco ' " 2 Jo (l + xx^/2)4 ^ ' ^ ^ 

In two interesting limits, energy loss behaves quite differently. At rj = (p = we have [11] 

de 2ashCp{Z,eBfe^ 

-dl = 3^^ ' (1'") 

^ = -0.37 a, h-^/^Cp iZgeBef/\ X»l- (12b) 

In the strong field limit energy loss is independent of the quark mass, whereas in the weak field 
case it decreases as Since x h, limit of x 1 corresponds to the classical energy loss. 

To apply this result to heavy-ion collisions we need to write down the invariant X ^ suitable 
kinematic variables. The geometry of a heavy-ion collision is depicted in Fig. 2. We can see that 
the vector of magnetic field B is orthogonal to the "reaction plane" , which is spanned by the impact 
parameter vector b and the collision axis (z-axes) . For a quark of momentum p we define the polar 
angle 6 with respect to the z-axis and azimuthal angle f with respect to the reaction plane. In 
this notation, B = By and p = PzZ + p± (x cos tp + ysmip), where p± = |p| sin Ri £ sin 6. Thus, 
{B X p)'^ = B^{p1 +p'^cos^ip). Conventionally, one expresses the longitudinal momentum and 
energy using the rapidity ry as e = m± coshr/ and Pz = m^sinhry, where = We have 

X^ = ^^pi(sinh2^ + cosV) (13) 
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FIG. 2: Geometry of a heavy-ion collision, p denotes the quark momentum. Two orthogonal planes are the 
"reaction plane" in which the initial heavy-ion momenta lie, and the mid-rapidity plane 9 = 7r/2, which is 
labeled as ?7 = 0. 

In Fig. 3 we show the calculation of the energy loss per unit length in a constant magnetic field 
using (11) and (13). We see that energy loss of a n quark with p± = 10 GeV is about 0.2 GeV/fm 
at RHIC and 1.3 GeV/fm at LHC. This corresponds to the loss of 10% and 65% of its initial energy 
after traveling 5 fm at RHIC and LHC respectfully. Therefore, energy loss due to the synchrotron 
radiation at LHC gives a phcnomcnologically important contribution to the total energy loss. 

Energy loss due to the synchrotron radiation has a very non-trivial azimuthal angle and rapidity 
dependence that comes from the corresponding dependence of the x-parameter (13). As can be 
seen in Fig. 3(c), energy loss has a minimum at (/? = 7r/2 that corresponds to quark's transverse 
momentum p± being parallel (or anti-parallel) to the field direction. It has a maximum at = 0, tt 
when pj_ is perpendicular to the field direction and thus lying in the reaction plane. It is obvious 
from (13) that at mid-rapidity 77 = the azimuthal angle dependence is much stronger pronounced 
than in the forward/backward direction. Let mc emphasize, that the energy loss (11) divided by 
m^, i.e. de/{dlm?) scales with x- In turn, x is a function of magnetic field, quark mass, rapidity 
and azimuthal angle. Therefore, all the features seen in Fig. 3 follow from this scaling behavior. 
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FIG. 3: Energy loss per unit length by quarks moving in constant external magnetic field as a function of (a) 
at RHIC at 77 = = 0, (b) -p^ at LHC at 77 = = 0, (c) azimuthal angle with respect to the reaction 
plane for prp = 5 GeV and 77 = 0, 1. (d) Fractional energy loss vs rapidity 77 at pt = 5 GeV and ip — 0, n/i. 

IV. POLARIZATION OF LIGHT QUARKS 



Synchrotron radiation leads to polarization of electrically charged fermions; this is known as 
the Sokolov-Ternov effect [9] . Unlike energy loss that we discussed so far, this is a purely quantum 
effect. It arises because the probability of the spin- flip transition depends on the orientation of the 
quark spin with respect to the direction of the magnetic field and on the sign of fermion's electric 
charge. The spin- flip probability per unit time reads [9] 

3 / 



w 



16" 



ZqC \v X B\ 



1 - I (C • vf - ^sign(e,) (C • ^) J , (14) 



where C is a unit axial vector that coincides with the direction of quark spin in its rest frame, 
V = p/e is the initial fermion velocity. 

The nature of this spin flip transition is transparent in the non-relativistic case, where it is 
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induced by the interaction [10] 



geZqh\ ^ ^ 



H = -i2.B = -(^-^] s-B, (15) 



It is seen, that negatively charged quarks and anti-quarks (e.g. u and d) tend to ahgn against the 
field, while the positively charged ones (e.g. u and d) align parallel to the field. 

Let Ti^d) be the number of fermions or anti-fermions with given momentum and spin direction 
parallel (anti-parallel) to the field produced in a given event. At initial time t = to the spin- 
asymmetry defined as 

A = (16) 

vanishes. Eq. (14) implies that at later times, a beam of non-polarized fermions develops a finite 
asymmetry given by [9] 

8\/3 



A 

15 

where 



(l-e-^), (17) 



r=— ^ i-Y ^^] (18) 

5V3masCF Ve/ \Zge\vxB\J 

is the characteristic time (do not confuse with (7)!) over which the maximal possible asymmetry is 
achieved. This time is extremely small on the scale of to- For example, it takes only r = 6 • 10^^ fm 
for a d quark of momentum p± = 1 GeV at r/ = 99 = at RHIC to achieve the maximal asymmetry 
of Ajn = = 92%. Therefore, quarks and anti-quarks are polarized almost instantaneously after 
being released from their wave functions. Note, that for quarks moving nearly parallel to the field 
direction, r can be large. However, the corresponding azimuthal angle and rapidity are always 
much smaller than their experimental resolution. 

As the quarks propagate through the hot medium their polarization is washed-out. This however 
does not seem to be an important effect because (i) the quark-gluon plasma is chirally symmetric 
phase of QCD and (ii) the spin-field interaction energy, estimated for, say, d quarks from (15) is 
of order 1.5 GeV, while the temperature in the plasma is of order of 0.2 GeV. Nevertheless, we do 
expect that the polarization of quarks is completely washed out in the fragmentation process at 
T <Tc because the chiral symmetry is broken in the QCD vacuum. 

The Sokolov-Ternov effect can be observed by studying the polarization of light charged leptons. 
The characteristic time r is still very small. Indeed, in (18) one replaces agCp by aem which - in 
the case of electron - is compensated by smaller mass. In my opinion, observation of such a lepton 
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polarization asymmetry would be the most definitive proof of existence of the strong magnetic field 
at early times after a heavy-ion collision regardless of its later time-dependence. 

V. CONCLUSIONS 

In this paper we argued that the magnetic field created by fast heavy ions can be considered 
as approximately constant due to high electric conductivity of the quark-gluon plasma. We then 
used some well-known QED results to estimate the energy loss suffered by fast quarks in external 
magnetic field in heavy-ion collisions. Wc found that the energy loss per unit length for a light 
quark with px = 15 GeV is about 0.27 GeV/fm at RHIC and 1.7 GeV/fm at LHC, which is 
comparable to the losses due to interaction with the plasma. Thus, the synchrotron radiation 
alone is be able to account for quenching of jets at LHC with as large as 20 GeV. Synchrotron 
radiation seems to be one of the missing pieces in the puzzle of the jet energy loss in heavy ion 
collisions. 

Wc also pointed out that quarks and Icptons arc expected to be strongly polarized in plasma in 
the direction parallel or anti-parallel to the magnetic field depending on the sign of their electric 
charge. While polarization of quarks is expected to be washed out during the fragmentation, 
that of leptons should survive to present a direct experimental evidence for the existence of the 
magnetic field. Finally, we would like to mention a possible deep connection between the "tunneling 
through the horizon" thermalization mechanism that we discussed before [25] and the polarization 
of fermions in external magnetic field that can be viewed as an Unruh effect [26-28] . 
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